HOC  FILE  COPY, ^  ADA091859 


ASS,  CATION  or  THIS  PAGI  IWim  Dalt  (alan«  JH  8® 

REPORT  DOCUMENTATION  PAGE  Ra*  1  P^Jl^EeVomplet^Ig^orm 

T^SWWniWlEir  ]i  OOVT  ACCESsioiTNa  1  RECIPIENT’*  CAT  AlOO  NUMBER 

c-t WE^fttrtranwwwT—  —  ...  .  . . . . s  tym-of  rcporj  »  period  co'  ereo 

i^y  C.iputer  Simulations  for  the  Hi  croscop^T?  .✓’’"X  /TT  Fl2al 
1  °f  Hot  Eros  I  on  .  , .  ^  VlTi  ■ 


Bruce  C ./Garrett  j  ^ 

t.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Battel  le  Columbus  Laboratories  /V 
Columbus,  OH  1*3201  ' 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ” 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 
IA  NOnTtorTnG  Aftrm-v  NAME  *  tnfumfK  f 

QMj 

I*.  DISTRIBUTION  STATEMENT  (ol  thit  Rwporl) 


cajH*ACT  or  grant  numbert*; 

T  DAAG29-77-C"00l*7 


to.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  »  WORK  UNIT  NUMBERS 


12.  RCPOBXOAT  E 


C / /  j  Sep  80  ) 

V~  *  -nr-  NUMBER  OF 


»  TU 

n$  Otllco)  15.  S  E  C  U R~l T Y  CLASS,  (ot  thtm  report) 

Unclassified 

1  sir  DECL ASSl FI CATI ON/ DOWN GR ADI NC- 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited, 


17.  DISTRIBUTION  STATEMENT  (ot  the  oboUoct  onlararf  In  Btock  20 ,  It  dltimrmnt  from  Report) 


IB.  SUPPLEMENTARY  NOTES 


The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


19.  KEY  WORDS  YConl/nu*  an  rcrara*  aide  tt  n#c •••ary  and  idontny  oy  oioc*  number) 


gases 

metals 

computerized  simulation 
e  ros  i  on 


surface  properties 
chemi cal  react i ons 


1 20  ABSTRACT  fCmmtamm  mm  »•* 


[  identity  by  block  tumibor) 


y  ^  systematic  technique  for  obtaining  physical  reasonable  models  of  gas-surface 
interaction  potentials  was  employed.  A  general  computer  model  has  been 
developed  to  provide  realistic,  interaction  potentials  for  chemical  reactions 
of  gas  species  on  a  variety  of  surface  structures.  In  order  to  better  under¬ 
stand  the  model  interaction  potentials  a  computer  program  was  developed  to 
help  probe  and  characterize  the  potential  energy  hypersurfaces.  This  enabled 
critical  configurations  such  as  equilibrium  geometries  and  saddle  points  to  be 


do 


EDITION  OF  •  MOV  (*  IS  OBSOLETE 


UNCLASSIFIED 


ri  »wirir*Tift»*  or  tu*«  »Acr  /■ 


ria«B 

Hr 


Unciacsl ilea 


SKCUNITV  CLASSIFICATION  OF  THIS  PAGEPThwi  0«<«  Bnfnd) 


20.  ABSTRACT  CONTINUED 


located;  the  program  also  provided  minimum  energy  pathways  for  reactions 
occurring  on  the  metal  surfaces.  Classical  trajectory  studies  were  carried 
out  for  models  of  direct  chemical  reaction  of  gas  phase  species  with  surfaces 
and  for  atomic  recombination  reactions  occurring  on  the  surfaces.  The  atomic 
recombination  studies  dealt  with  hydrogen  and  nitrogen  atom  recombination  on 
Fe  and  Ta  surfaces.  These  studies  showed  that  the  heat  inputed  into  the  surfac 
from  the  atomic  recombination  reaction  can  be  a  significant  fraction  of  the 
total  energy  of  the  heat  release  of  the  process. 


'Acce'r>r,'°n  >r  - 
i  »r*IC 


Avoi'i**-5'  _  •  /0f 


Ift: 


FINAL  REPORT 


on 

Contract  #  DAAG29-77-C-0047 

COMPUTER  SIMULATIONS  FOR  THE  MICROSCOPIC 
MECHANISM  OF  HOT  GAS  EROSION 

to 

U.  S.  Army  Research  Office 

September  30,  1980 


Prepared  by 


Bruce  C.  Garrett 


BATTELLE 

Columbus  Laboratories 
505  King  Avenue 
Columbus,  Ohio  43201 


TABLE  OF  CONTENTS 


Page 


I.  PROJECT  GOALS . 1 

II.  PROJECT  ACCOMPLISHMENTS . 3 

A.  Summary . 3 

B.  Details . 3 

C.  Publications  Resulting  from  this  Project . 5 

0.  Scientific  Personnel.  Supported  on  this  Project . 5 

III.  TECHNICAL  DISCUSSION . 6 

A.  Modeling  of  the  Potential  Energy  Surface . 6 

B.  Classical  Trajectory  Methods . 12 

REFERENCES . 18 


LIST  OF  TABLES 


Pag  a 

TABLE  I.  ATOM-ATOM  MORSE  PARAMETERS  AND  LATTICE 

PARAMETERS  FOR  THREE  METAL  SURFACES  .  19 

TABLE  II.  LATTICE  POTENTIAL  PARAMETERS  FOR  FIVE 

SURFACES  OF  Fe  . 19 

TABLE  III.  LATTICE  POTENTIAL  PARAMETERS  FOR  FIVE 

SURFACES  OF  7  20 

G 

TABLE  IV.  LATTICE  POTENTIAL  PARAMETER  FOR  FIVE 

SURFACES  OF  Cr  . ' . .  .  10 

TABLE  V.  ATOM-ATOM  MORSE  PARAMETERS  FOR 

GAS  PHASE  DIATOMICS . 21 

TABLE  VI.  ATOM- ATOM  MORSE  PARAMETERS  FOR  GAS-ATOM 

INTERACTIONS  WITH  SEMI-INFINITE  SURFACES . 21 

TABLE  VII.  HYDROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Fe . 22 

TABLE  VIII.  HYDROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Ta  . 22 

TABLE  IX.  HYDROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Cr  . 23 

TABLE  X.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Fe  . 24 

TABLE  XI.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Ta  . 25 

TABLE  XII.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE 

SURFACES  OF  Cr  . 26 

TABLE  XIII.  Fe-Fe  MORSE  POTENTIAL  PARAMETERS  FOR 

A  FINITE  CLUSTER  OF  Fe  ATOMS . 27 

TABLE  XIV.  ATOM-ATOM  MORSE  PARAMETERS  FOR  HYDROGEN 

ADSORPTION  ON  FINITE  Fe  CLUSTERS . 28 

TABLE  XV.  ATOM-ATOM  MORSE  PARAMETERS  FOR  HYDROGEN 

ADSORPTION  ON  FINITE  Ta  CLUSTERS . 2G 


LIST  OF  TABLES 
(Continued) 


Pages 


TABLE  XVI.  ATOM-ATOM  MORSE  PARAMETER  FOR  HYDROGEN 

ADSORPTION  ON  FINITE  Cr  CLUSTERS . 29 

TABLE  XVII.  ATOM-ATOM  MORSE  PARAMETERS  FOR  NITROGEN 

ADSORPTION  ON  FINITE  Fe  CLUSTERS  .  30 

TABLE  XVIII.  ATOM-ATOM  MORSE  PARAMETERS  FOR  NITROGEN 

ADSORPTION  ON  FINITE  Ta  CLUSTERS . 31 

TABLE  XIX.  ATOM-ATOM  MORSE  PARAMETERS  IN  NITROGEN 

ADSORPTION  ON  FINITE  Cr  CLUSTERS . 32 

TABLE  XX.  ENERGETICS  OF  ATOMIC  RECOMBINATION  ON 

METAL  SURFACES . 33 

TABLE  XXI.  TRAJECTORY  RESULTS  FOR  ATOM  RECOMBINATION  ON 

METAL  SURFACES  BY  A  RADEAL  MECHANISM . 34 

TABLE  XXII.  RESULTS  OF  TRAJECTORY  CALCULATIONS  FOR  NITROGEN 
ATOM  RECOMBINATION  ON  Fe  SURFACES  FOR  NON-NORMAL 
INCIDENCE  OF  THE  GAS  ATOM . 36 


rranfr 


« *! 


I. 


PROJECT  GOALS 


An  understanding  of  the  mechanisms  of  the  erosion  of  metal  surfaces 
in  hot  flowing  media  is  important  in  assessing  methods  of  lessening  the  effects 
of  wear  and  erosion  upon  the  performance  of  large  caliber  gun  tubes.  The 
overall  mechanisms  of  erosion  are  extremely  complex  including  hydrodynamic 
flow,  chemical  reactions,  and  direct  erosion  by  macroscopic  particles. 

Recent  experimental  studies^  indicate  that  chemical  processes  may  be 
important  in  these  mechanisms.  In  this  report  we  present  the  results  of 
theoretical  studies  of  the  microscopic  mechanism  of  chemical  erosion  which 
are  aimed  at  assessing  the  importance  of  such  processes  under  the  physical 
conditions  found  in  gun  tubes. 

We  consider  two  broad  possibilities  of  chemical  erosion:  (1)  chemical 

heating  of  the  surface  causing  local  melting  and  subsequent  removal  of  the 

material  by  the  flowing  gas  stream,  and  (2)  direct  chemical  erosion  in  which 

gaseous  species  attack  surface  atoms  forming  volatile  chemical  species 

leading  to  removal  of  the  surface  atoms.  By  chemical  heating  of  the  surface 

we  mean  heat  input  into  the  surface  from  chemical  processes  occurring  near 

the  surface.  We  therefore  exclude  simple  convective  heating  from  our 

consideration  although  it  is  undoubtably  important  in  the  erosion  process. 

At  least  two  distinct  mechanisms  of  chemical  heating  can  be 

discerned:  (1)  heating  from  chemisorption  and  (2)  heating  from  recombination 

reactions.  In  the  process  of  forming  bonds  between  gaseous  species  and 

the  surface,  energy  can  be  released  into  the  surface.  For  example,  when  0£ 

3  2 

chemisorbs  on  steel  it  releases  4.7  x  10  kO/kg-metal .  When  two  reactive 

species  such  as  free  radials  recombine  to  form  a  single  molecule  heat  can 

3 

be  released  into  the  surface. 
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Direct  chemical  reaction  of  gaseous  species  with  the  surface  to 
form  a  gas-phase  species  containing  a  metal  atom  is  similar  to  the  process 
of  chemisorption  heating  and  the  subsequent  wipe-off  of  the  metal  atom  by 
the  flowing  gas.  In  both' cases  the  energy  liberated  from  the  formation 
of  a  gas-metal  bond  is  used  to  severe. the  bending  between  the  metal  atom 
and  the  remainder  of  the  surface.  The  major  difference  is  that  in 
direct  chemical  erosion  the  metal  atom  leaves  the  surface  in  a  bound 
complex  with  the  gas  atom.  In  direct  chemical  erosion  a  surface  atom 
is  removed  in  a  single  collision  event,  whereas,  in  the  heat  and  wipe-off 
mechanism  many  collisions  can  heat  the  surface  leading  to  the  removal  of 
surface  atoms. 

In  summary,  the  goals  of  this  project  were  the  investigation  of 
chemical  mechanisms  of  hot  gas  erosion  of  metal  surfaces  using  theoretical 
techniques.  In  particular,  the  processes  of  interest  are:  (1)  chemical  heatinq 
of  metal  surfaces  and  (2)  direct  chemical  erosion  of  metal  surfaces.  These 
studies  were  carried  out  for  a  variety  of  surface  structures  including 
perfect  crystal  surfaces  and  also  surfaces  with  defects  such  as  steps, 
kinks,  and  adatoms.  We  investigated  the  importance  of  the  surface  defects 
in  the  processes  mentioned  above, 
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II. _ PROJECT  ACCOMPLISHMENTS 

A.  Summary. 

A  systematic  technique  for  obtaining  physical  reasonable  models  of 
gas-surface  interaction  potentials  was  employed.  A  general  computer  model 
has  been  developed  to  provide  realistic,  interaction  potentials  for  chemical 
reactions  of  gas  species  on  a  variety  of  surface  structures. 

In  order  to  better  understand  the  model  interaction  potentials  a 
computer  program  was  developed  to  help  probe  and  characterize  the  potential 
energy  hypersurfaces.  This  enabled  critical  configurations  such  as  equilibrium 
geometries  and  saddle  points  to  be  located;  the  program  also  provided 
minimum  energy  pathways  for  reactions  occurring  on  the  metal  surfaces. 

Classical  trajectory  studies  were  carried  out  for  models  of  direct 
chemical  reaction  of  gas  phase  species  with  surfaces  and  for  atomic  recombina¬ 
tion  reactions  occurring  on  the  surfaces.  The  atomic  recombination  studies 
dealt  with  hydrogen  and  nitrogen  atom  recombination  on  Fe  and  Ta  surfaces. 

These  studies  showed  that  the  heat  inputed  into  the  surface  from  the  atomic 
recombination  reaction  can  be  a  significant  fraction  of  the  total  energy  of 
the  heat  release  of  the  process. 

B.  Details. 

In  this  work  we  studied  interactions  of  atomic  and  diatomic  gas 
species  with  surfaces.  Of  particular  interest  was  the  transfer  of  heat  from 
reacting  gas  species  into  the  surface.  Therefore,  the  model  had  to  allow  for 
motion  of  the  surface  atoms  as  well  as  the  gas-phase  species.  Such  a  model 


i 


S?  -LW .  M1^*?^** 


4 


4  5 

was  developed  and  has  been  previously  described.  *  This  model  is  semiempirical 
and  the  parameters  of  the  model  were  chosen  to  best  reproduce  experimental 
data  when  it  was  available.  An  important  feature  of  ;'ivs  model  is  its  ability 
to  model  surfaces  with  steps,  kinks,  adatoms,  etc.  .  her  details  of  the 
surface  fitting  for  parti  cl ar  systems  are  presented  in  the  technical  section. 

Direct  chemical  erosion  of  metal  surfaces  can  occur  if  sufficient 
energy  from  chemisorption  of  the  gaseous  species  is  transferred  into  breaking 
the  bond  between  a  metal  atom  and  the  rest  of  the  surface.  Surface  atoms  at 
step-sites,  kink-sites,  or  surface  adatoms  are  bound  more  weakly  to  the 
surface  than  surface  atoms  with  their  full  compliment  of  nearest  neighbors. 
Therefore,  one  might  expect  that  direct  chemical  erosion  of  surface  atoms  at 
step-  or  kink-sites  may  be  more  favorable  than  for  surface  atoms  at  perfect 
crystal  sites.  However,  we  have  run  several  thousand  trajectories  on  surfaces 
with  various  defects  -  i.e.,  step,  kink,  etc.  -  and  have  found  no  evidence 
of  direct  chemical  erosion. 

Studies  of  recombination  heating  of  the  solid  by  a  Langmuir- 
Hinshelwood  mechanism  have  been  previously  reported.6  We  report  here  work 
in  which  we  have  studied  in  detail  recombination  occurring  by  a  Rideal 
mechanism.  We  have  also  examined  the  effects  of  steps  and  kinks  upon  the 
heating  of  the  surface.  A  surprising  outcome  of  this  study  was  the  relative 
insensitivity  of  the  heating  of  the  surface  upon  parameters  such  as  tempera¬ 
ture  of  the  solid  and  incident  translational  energy  of  the  impringing  atom. 
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TECHNICAL  DISCUSSION 


A.  Modeling  of  the  Potent* a 1  Energy  Surface. 


The  model  potential  used  for  the  interaction  of  diatomic  gas 

species  with  solid  surfaces  is  described  in  Appendix  A.  The  potential  VLEps 

is  based  upon  the  semiempirical  London-Eyring-Polanyi-Sato  (LEPS)  description 

of  gas-solid  bonding/  The  functional  form  is  given  in  Equations  5-9  of 

Appendix  A.  The  parameters  which  must  be  specified  are  the  dissociation 

energy  Dpg,  the  Morse  range  parameters  apg,  the  equilibrium  geometry  Rpg 

and  the  Sato  parameter  Apg  for  the  two  body  interactions  between  atoms  P  and 

Q.  P  (and  Q)  can  be  A  or  8,  one  of  the  gas  atoms,  or  one  of  the  surface  atoms. 

This  model  is  capable  of  describing  the  interaction  of  a  diatomic  mc"!ecu1e 

with  a  surface  of  moving  atoms  -  the  instananeous  locations  of  the  surface 

atoms  are  used  in  evaluating  the  potential  expression  -  however,  the  total 

potential  energy  is  the  sum  of  V^ps  and  the  potential  energy  obtained 

from  displacements  of  the  lattice  from  its  equilibrium  geometry.  The  lattice 

potential  was  described  by  a  modified  Einstein  model  in  which  the  potential 

was  given  in  terms  of  displacements  of  the  atoms  from  their  individual 

equilibrium  positions  in  a  rigid  lattice.  Because  we  were  interested  in 

describing  situations  in  which  large  amounts  of  energy  can  be  transferred 

to  the  surface  we  used  an  anharmonic  potential,  the  Morse  potential, 

8 

rather  than  the  usual  harmonic  Einstein  model.  In  the  following  we 
describe  how  the  parameters  of  this  model  were  determined.  In  all 
cases  the  parameters  were  calculated  so  that  the  models  reproduced 
experimental  data  when  it  was  available. 
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Consider  first  the  parameters  for  the  lattice  potential 


VLAT  =  DLAT  ^ 


.e'^ATd]2  - 


LAT 


(1) 


where  d  is  the  displacement  from  the  equilibrium  lattice  site  location. 
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Girifalco  and  Weizer  used  the  sum  of  two- body  Morse  potential  between  all 
surface  atoms  to  describe  metal  bonding.  Using  the  results  of  their  work 
we  obtained  the  parameters  of  the  modified  Einstein  oscillators.  From 
Girifalco  and  Weizer^  work  the  total  potential  $  felt  by  one  "test  atom" 
in  the  lattice  was  obtained  by  summing  over  the  two-body  interactions  between 
the  "test  atom"  and  all  other  atoms  in  the  lattice.  The  lattice  was  assumed 
to  be  semi -infinite,  that  is  the  sums  were  over  lattice  positions  with  x  and  y 
ranging  from  -®  to  ®  and  z  from  0  to  -®.  In  practice  the  sums  were  done 
over  a  sufficient  number  of  lattice  sites  to  insure  convergence.  The  "test 
atom"  was  taken  to  be  in  the  top  layer  of  the  lattice;  all  other  atoms  in 
the  lower  layers  were  fixed  at  their  equilibrium  lattice  sites.  The 
lattice  constant  (a)  was  taken  from  Girifalco  and  Weizer  when  available  or 
from  Pea~son^and  we  considered  only  100  faces  of  the  lattices.  The  locations 
of  the  atoms  in  the  top  layer  were  relaxed  from  their  lattice  sites  so  that 
the  forces  on  the  test  atom  vanished.  A  more  physical  model  would  also 
include  relaxation  of  the  lower  layers,  however,  since  this  effect  was  small 
for  the  top  layer  and  decreased  for  the  lattice  layers  below  the  surfaces  we 
neglected  it  for  all  but  the  top  layer.  This  determined  the  equilibrium 
geometry  and  therefore  the  value  of  the  potential  at  this  location  was  the 
negative  of  the  energy  necessary  to  remove  the  atom  from  the  surface. 


(This  assumes  the  zero  of  energy  for  all  the  two-body  potentials  was  at 


infinite  separation  of  the  atoms.)  This  defined  D^.  and  the  range  parameters 
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evaluated  at  the  equilibrium  geometry.  The  z  direction  was  taken  to  be 
perpendicular  to  the  surface,  the ref ore, c,  was  the  range  parameter  for  motion 

LM  1 

normal  to  the  surface.  In  general,  the  range  parameters  for  motion  in 
the  surface  plane  will  be  different  than  for  out  of  the  plane,  but  we  did 
not  take  this  into  account  in  our  model. 

Using  the  sum  of  pair  potentials  we  assessed  the  effect  of  surface 
defects  (steps,  kinks,  etc.)  upon  the  bonding  of  the  metal  atoms  to  the 
surface.  As  surface  atoms  were  removed,  the  "test  atom"  became  less 
tightly  bound  decreasing  D^y,  and  also  the  equilibrium  geometry  and 
range  parameters  for  the  test  atom  changed.  We  considered  five  types  of 
surfaces:  (1)  a  regular  crystal  surface  (no  defect),  (2)  a 
simple  step  in  which  all  the  top  layer  atoms  in  the  left  half  of  the  x-y 
plane  had  been  removed,  (3)  a  kink  in  which  all  the  top  layer  atoms  in 
all  but  the  upper  right  quadrant  of  the  x-y  plane  had  been  removed, 

(4)  an  adatom  next  to  a  simple  step,  and  (5)  a  lone  adatom  (removal  of 
all  the  top  surface  atoms  except  for  the  "test  atom.") 

In  Table  I  the  atom-atom  Morse  parameters  from  Girifalco 
and  Weizer's  work  along  with  the  lattice  parameters  for  Fe,  Ta,  and  Cr 
surfaces  are  listed.  In  Tables  1 1 - IV  we  present  our  calculated  equilibrium 
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geometries  and  lattice  potential  parameters  for  these  five  surface  defects 
on  Fe,  Ta,  and  Cr  surfaces.  As  the  number  of  nearest  neighbors  decreas  d 
the  bonding  of  the  "test  atom"  also  decreased.  The  sublimatio  energy 
at  zero  temperature  and  pressure  is  a  weighted  average  of  the  dissociation 
energies  for  the  difference  surface  sites  of  the  lattice,  therefore  we 
expected  our  values  of  for  the  difference  surface  defects  to  bracket  the 
experimental  sublimation  energy.  For  Fe  and  Cr  the  sublimation  energies 
are  97  and  90  kcal/mol,  respectively!1  and  for  Ta  it  is  estimated  from 
the  experimental  values  of  Mo  and  W  to  be  between  160  and  200  kcal/mol. 

We  next  consider  the  manner  in  which  the  parameters  of  the  gas- 
surface  interaction  potential  VLEps  were  determined.  In  the  limit  that 
the  gas-phase  atom  A  and  B  are  well  above  the  surface, VEEps  reduces  to  a 
Morse  function  with  dissociation  energy  DAB>  range  parameter  aAB>  and 
equilibrium  geometry  RAg.  These  values  were  fit  from  experimental  values 
for  Hg  and  and  for  FeH  we  used  a  theoretical  estimate. ^>14 

'Values  of  these  parameters  are  listed  in  Table  V. 

In  the  limit  that  one  of  the  gas-phase  atoms  A  is  far  removed 
from  the  surface  and  the  other  B  is  on  the  surface,  the  potential  equals 
the  adsorption  energy  of  atom  B  on  the  surface.  This  was  obtained  by 
summing  the  pair  potential  between  the  adsorbing  atom  and  all  of  the 
surface  atom.  The  pair  potential  is  a  Morse  function  with  dissociation 
energy  Dg^,  range  parameters  a^,  and  equilibrium  geometry  Rgg,  where  S 
denotes  a  surface  atom.  Olander'  has  treated  the  adsorption  of  hydrogen  atoms 
on  Fe  and  Ta  surfaces  of  infinite  lattices  without  surface  defects. 
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The  gas-surface  Morse  parameters  for  H  on  Cr  and  for  N  on  all  three  metal 
surfaces  have  been  estimated.^  The  atom-atom  Morse  parameters  for  these 
systems  are  presented  in  Table  VI.  The  location  of  the  adsorbing  atom 
was  optimized  such  that  the  forces  on  it  vanish.  The  potential  at  this 
location  was  the  negative  of  the  adsorption  energy,  and  we  also  calculated 
the  frequency  of  the  motion  of  the  adsorbed  atom  normal  to  the  surface. 
These  values  are  listed  in  Tables  VT I -XII  for  H  on  N  adsorbing  on  three  of 
the  surface  types  {no  defect,  step,  and  kink)  for  Fe,  Ta,  and  Cr. 

In  our  dynamical  calculations  it  was  impractical  to  use  semi- 

infinite  surfaces,  therefore,  it  was  necessary  to  consider  small  clusters 
of  metal  atoms,  30-46  atoms.  In  the  sum  of  pairs  model,  if  the  number 
of  surface  atoms  was  decreased  the  bonding  of  the  adsorbed  atom  also 
decreased.  Therefore,  we  found  it  necessary  to  reoptimize  the  atom-atom 
interaction  parameters  so  that  the  adsorption  energy  and  vibrational 
frequency  of  the  adsorbed  atom  would  be  the  same  as  on  the  semi -Infinite 
surfaces.  This  was  done  by  finding  the  values  of  Dgs,  ags,  and  Rgs  that 
would  simultaneously  give  the  correct  equilibrium  height  above  the  surface, 
adsorption  energy, and  vibrational  frequency  for  the  adsorbing  atom.  These 
parameters  depend  upon  the  size  of  the  cluster.  We  found  that  different 
cluster  sizes  were  needed  to  properly  describe  the  different  dynamical 
events  -  direct  chemical  reaction  and  recombination  heating  -  and  we 
describe  both  below. 

First  consider  the  direct  chemical  reaction  of  a  gas  atom  A 
with  a  surface  atom  to  form  a  gas-phase  species  in  which  A  is  bonded  to 
the  metal  atom.  We  treated  the  metal  atom  that  was  attacked  as  the 
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species  B  adsorbed  to  the  surface.  The  parameters  for  the  two-body 
potential  between  B  and  the  other  surface  atoms  was  just  those  given 
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by  Girifalco  and  Weizer  for  the  case  of  an  semi-infinite  lattice. 

In  our  dynamical  calculations  we  used  a  cluscer  of  45  metal  atoms  plus 
the  atom  B  to  represent  the  surfaces  with  no  defects.  For  a  surface  with 
no  defects,  twenty  atoms  were  used  in  the  top  layer  (a  5  x  5  grid  with 
atoms  removed  from  the  corner  and  the  central  atom  replaced  by  atom  B), 

12  in  the  second  layer  (a  4  x  4  grid  with  atoms  removed  from  the  corners), 

9  in  the  third  layer  (a  3  x  3  grid),  and  4  is  the  fourth  layer  (a  2  x  2  grid). 

The  reoptimized  atom-atom  Morse  parameter  for  this  finite  lattice  are 
presented  in  Table  XIII  for  the  five  surface  defects  described  above. 

These  parameters  should  be  compared  with  those  in  Table  I  for  an  estimate 
of  the  effect  of  finite  cluster  size  upon  the  Morse  parameters. 

In  modeling  the  direct  chemical  erosion  process  we  used  the  Morse 
parameters  for  FeH  listed  in  Table  V  and  for  H  on  the  Fe  surface  v:e  used 
the  parameters  given  in  Table  XIV.  The  specification  or  the  interaction 
potential  VL£PS  is  completed  by  selecting  the  Sato  parameters.  The  Sato 
parameters  were  varied  to  obtain  the  best  qualitatively  correct  description  of  the 
adsorption  of  H  on  the  Fe  surface.  We  found  that  for  all  a ' s  equal  to 
zero  the  adsorption  energy  on  defect  free  Fe  was  50  kcal/mol  with  an 
equilibrium  height  of  H  above  the  surface  of  3.6  aQ.  These  numbers  compare 
well  with  the  values  in  Table  VII, 

For  describing  recombination  heating  we  used  a  cluster  size  of 
30  atoms  for  the  lattice  with  no  surface  defect.  This  cluster  had  16 
atoms  in  the  top  layer  (a  4  x  4  grid),  9  atoms  in  the  second  layer 
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(a  3  x  2  grid),  4  atoms  in  the  third  layer  (a  2  x  2  grid)  and  finally  a 
single  atom  in  the  fourth  layer.  The  reoptimized  atom-atom  Morse  parameters 
for  H  and  N  adsorbtion  on  Fe,  Ta,  and  Cr  are  summarized  in  Tables  XIV-XIX 
for  three  types  of  surfaces  -  no  defect,  step,  and  kink.  These  should 
be  compared  with  those  in  Table  VI  to  see  the  effect  of  finite  cluster 
size  upon  the  potential  parameters  describing  adsorbtion. 

The  Morse  parameters  used  in  V^pS  can  be  found  in  Tables  V 
and  XIV-XIX.  The  specification  of  the  interaction  potential  is  completed 
by  selecting  values  for  the  Sato  parameters.  Because  we  had  no  guide  for 
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their  selection  we  used  the  values  zero  as  in  previous  applications.  We 
performed  a  detailed  study  of  the  interaction  potential  used  for  the 
Rideal  recombination  of  H2  on  Fe.  The  adsorbed  H  was  placed  at  a  4-fold 
equilibrium  site  3.13  aQ  above  the  top  layer.  The  second  H  atom  was 
brought  down  directly  over  the  adsorbed  H  and  the  adsorbed  H  was  allowed 
to  relax  to  its  equilibrium  position  and  a  local  minimum  in  the  potential 
was  found  when  the  second  H  was  4.76  aQ  above  the  surface  and  the 
adsorbed  H  was  2.45  aQ  above  the  surface.  This  local  minimum  was  higher 
in  energy  by  only  0.4  kcal/mol  than  desorbed  H2  far  removed  from  the 
surface.  However,  there  is  a  6.3  kcal/mol  barrier  to  this  desorption 
process.  The  barrier  to  this  process  had  its  maximum  when  the  desorbing 
atoms  were  6.63  and  5.14  kcal/mol  above  the  surface. 

B.  Classical  Trajectory  Methods. 

All  of  the  trajectory  calculations  performed  here  were  for 
the  general  event  in  which  a  single  atom  collided  with  an  adsorbed  atom 
and  the  desired  outcome  was  a  reactive  event  in  which  the  adsorbed  atom 
and  incident  gas  atom  form  a  new  bond  and  leave  the  surface.  A 
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sufficiently  large  cluster  of  metal  atoms  was  used  to  give  a  physically 

correct  description  of  the  interaction  potential.  Only  a  small  number, 

generally  4-8,  of  the  surface  atoms  were  allowed  to  move.  The  others 

were  fixed  at  their  lattice  site  locations.  A  unique  trajectory  was 

specified  by  selection  of  the  'Initial  position  and  momentum  vectors 

of  the  moving  surface  atoms,  the  adsorbed  atom,  and  the  incident  atom. 

It  was  assumed  that  the  surface  and  adsorbed  atom  are  at  thermal  equilibrium 

characterized  by  a  temperature  T.  The  momenta  for  these  atoms  were 

selected  from  a  Boltzmann  distribution  by  importance  sampling.^  The 

positions  of  these  atoms  were  also  selected  from  a  Boltzmann  distribution 

18 

but  using  Metropolis  sampling.  The  incident  atom  was  started  at  an 
arbitrary  and  large  distance  above  the  surface,  generally  10  aQ.  The 
incident  angle  was  specified  and  the  x,y  coordinates  were  chosen  so  that 
an  undeflected  straight  line  trajectory  would  hit  the  surface  at  xQ,  y  . 

The  coordinates  xQ,  yQ  were  chosen  randomly  on  a  square  of  length  equal  to 
the  lattice  parameter  .  The  momentum  vector  was  chosen  along  the  incident 
angle  with  magnitude  equal  to  the  specified  initial  translation  energy 
of  the  atom  plus  any  potential  attraction  it  may  feel  from  the  surface. 

This  potential  term  was  evaluated  by  subtracting  the  value  of  the  total 
potential  with  the  incident  atom  at  its  initial  location  from  the  total 
potential  with  the  incident  atom  moved  infinitely  far  from  the  surface. 

Direct  chemical  erosion  was  modeled  only  for  H  incident  upon  Fe 
surfaces.  In  the  gas-phase  Fe-H  is  bound  by  about  60  kcal/mol  and  Fe 
is  bound  to  the  metal  surface  by  74-128  kcal/mol  depending  upon  whether 
the  Fe  atom  is  in  a  smooth  face  of  the  lattice  or  at  a  surface  defect  site. 
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Therefore  the  process  of  H  pulling  a  Fe  atom  off  the  surface  is  uphill 
in  energy  by  14-68  kcal/mol.  The  most  energetically  favorable  situation 
is  for  an  H  atom  to  attack  an  adatom  on  the  surface.  Even  in  this 
situation  the  H  atom  must  have  about  14  kcal/mol  (0.6  eV)  of  translational 
energy  to  surmount  the  barrier  to  form  FeH(g).  At  these  relatively  high 
translational  energies  the  light  hydrogen  atom  is  moving  very  rapidly 
making  the  prospects  of  it  interacting  strongly  with  a  heavy  Fe  atom  very 
low.  From  this  kinematic  argument  it  is  not  surprising  that  no  evidence 
was  seen  for  direct  chemical  erosion  with  hydrogen.  To  further  study 
this  problem  we  extended  our  model  calculations  to  include  collision  of 
heavier  incident  atoms  with  more  favorable  mass  ratios  with  the  Fe. 

The  potential  was  still  modeled  using  the  hydrogen  parameters,  and  only 
the  mass  of  the  incident  atom  was  changed  so  that  we  could  assess  whether 
the  nonreactivity  of  the  Fe  surface  was  due  to  this  kinematic  effect 
or  a  more  subtle  effect  from  a  combination  of  the  potential  surface 
and  structure  of  the  surface.  When  we  increased  the  mass  of  the  incident 
atom  we  still  saw  no  evidence  of  direct  chemical  erosion  and  therefore 
we  believe  that  this  mechanism  will  contribute  negligibly  to  the  hot-gas 
erosiGn  of  Fe  surfaces.  Tantalum  atoms  are  even  more  strongly  bound 
to  the  surface  than  are  Fe  atoms  and  direct  chemical  reaction  with  Ta 
atoms  in  the  surface  is  energetically  less  favorable  than  with  Fe  surfaces. 
The  surface  structure,  sublimation  energy,  and  adsorption  of  hydrogen 
is  very  similar  for  Fe  and  Cr,  and,  therefore,  the  same  general  conclusion 
will  hold  for  this  surface  also. 
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We  have  performed  detailed  studies  of  recombination  heating 
by  a  Rideal  mechanism  for  a  subset  of  the  gas-surface  interactions  we 
have  modeled.  Hydrogen-atom  recombination  has  been  studied  on  defect 
free  a  Fe  surface.  Nitrogen-atom  recombination  on  Fe  surfaces  has 
been  studied  for  surfaces  with  no  defects  and  surfaces  with  steps  and 
kinks  for  the  two  different  sets  of  N-Fe  Morse  parameters.  Nitrogen-atom 
recombination  on  Ta  surfaces  has  been  studied  for  surfaces  with  no 
defects  and  with  steps.  The  atom  recombination  reactions  of  all  of  the 
systems  are  exoergetic,  that  is  they  release  energy;  the  net  exoergicity 
is  the  difference  between  the  energy  of  the  new  gas-phase  diatomic  bond 
that  is  formed  and  the  adsorption  energy  of  the  atom.  The  values  of  the 
exoergici ties  for  the  various  systems  are  given  in  Table  XX.  We  were 
interested  in  how  much  of  the  exoergicity  can  be  deposited  into  the  surface. 
Because  the  nitrogen  molecule  forms  a  much  stronger  bond  than  the  hydrogen 
molecule,  N-atom  recombination  is  a  better  choice  for  studying  the  surface 
heating.  In  fact,  for  hydrogen  recombination  on  Ta  surfaces  the  exoergicity 
is  only  6-10  kcal/mol  depending  upon  the  surface  structure  { i . e . ,  step 
or  kinks,  etc.).  This  is  the  reason  that  most  of  our  detailed  studies 
are  on  N-atom  recombination. 

Although  the  recombination  process  is  energetically  favorable 
for  all  of  these  systems,  for  some  of  the  systems  there  is  a  more  stable 
intermediate  complex  -  the  adsorption  of  two  atoms  at  different  equilibrium 
sites  on  the  surfaces.  This  energy  will  be  close  to  the  sum  of  adsorption 
for  the  two  atom  independent  but  slightly  different  because  of  their 
interaction  which  is  normally  repulsive  in  our  models.  In  all  our 
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models  the  diatomic  dissociated  upon  adsorption.  We  found  that  for 
on  Fe  and  N2  on  Fe  using  Morse  parameter  set  (2),  that  two  adsorbed 
atoms  were  the  most  favorable  geometry.  For  these  systems  the  reaction 
probabilities  are  small  on  a  short  time  scale,  however,  given  a 
sufficiently  long  time  the  probability  increases  that  the  two  atoms 
will  recombine  and  leave  the  surface.  However,  this  then  becomes  a 
Langmiur-Hinshelwood  type  mechanism  which  we  exclude  from  our  considera¬ 
tion  here.  Work  on  this  mechanism  has  been  treated  earlier.^ 

A  summary  of  our  trajectory  results  are  presented  in  Tables 
XXI  and  XXII.  As  noted  above  the  two  systems  H2  on  Fe  and  N2  on  Fe  with 
Morse  parameter  set  2  show  small  probabilities  for  a  direct  recombination 
reaction,  the  primary  event  to  occur  is  adsorption  of  the  incident  atom. 
The  estimate  of  energy  transfer  for  these  systems  is  crude  because  of 
the  small  number  of  reactive  events,  nonetheless,  the  estimates  are  that 
between  10%  and  50%  of  the  exoergicity  is  deposited  into  the  surface. 

For  the  other  systems  the  statistics  are  much  better.  For  N2  on  Fe 
using  the  Morse  parameter  set  1  and  for  normal  incidence  upon  the  surface, 
the  fraction  of  exoergicity  deposited  into  the  surface  is  about  24-33%, 
Independent  of  the  incident  energy,  surface  temperature,  or  surface 
structure.  For  N2  on  Ta  using  Morse  parameter  set  1  the  fraction  of 
exoergicity  deposited  into  the  surface  in  the  range  5-11%  and  again  it 
is  independent  of  incident  energy,  surface  temperature,  or  surface 
structure.  The  trajectory  results  for  gas  atoms  impringing  on  the 
surface  at  an  angle,  see  Table  XXII,  show  a  slight  enhancement  of  the 
fraction  of  energy  transferred  to  the  surface  by  the  results  are  similar 
to  those  from  the  calculation  with  normal  incidence. 
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The  major  conclusion  to  be  drawn  from  these  studies  are  that 
(1)  atom- recombination  heating  of  the  surface  can  be  quite  important— up 
to  40%  of  the  bond  energy  of  was  found  to  be  transferred  to  the  Fe 
surface  and  (2)  energy  transfer  to  the  surface  is  more  facile  for  Fe 
than  for  Ta  surfaces.  For  systems  in  which  the  probability  of  the  Rideal 
mechanism  is  low  recombination  can  still  occur  by  a  Langmiun-Hinshelwood 
mechanism.  For  this  type  of  system  we  expect  heating  of  the  surface  to 
be  less  important  since  a  great  deal  of  the  new  bond  energy  is  used  up 
in  desorption  from  the  surface.  Ta  is  a  heavier  atom  ard  is  more  cohesive 
in  a  lattice  than  Fe.  The  heavy  mass  and  tight  bonding  of  Ta  tend  to 
make  it  'less  perturbed  by  the  light  gas  atom  and  thus  it  adsorbs  less  energy. 

A  surprising  outcome  is  that  the  surface  defects  had  little  effect  upon  either 
the  reaction  probability  PR  or  the  energy  transfer.  However,  this  can  be 
viewed  as  an  encouraging  results  since  the  simplier  models  of  perfect 
crystal  surfaces  with  no  defects  give  the  right  qualitative  behavior. 
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TABLE  I.  ATOM-ATOM  MORSE  PARAMETERS  AND  LATTICE 
PARAMETERS  FOR  THREE  METAL  SURFACES 


Metal 

D^kcal/mol) 

aM^a0  ^ 

rM(a0) 

Aa„) 

Fe 

9.63 

0.735 

5.38 

5.42 

Tae 

21.50 

0.731 

5.97 

6.25 

Cr 

10.18 

0.832 

5.20 

5.44 

Morse  dissociation  energy, 
k  Morse  range  parameter. 
c  Morse  equilibrium  geometry. 
d  Lattice  parameter. 

e  These  values  were  estimated  from  the  results  for 
Mq  and  W  of  Reference  9  (see  Reference  15  ). 


TABLE  II.  LATTICE  POTENTIAL  PARAMETERS  FOR  FIVE 
SURFACES  OF  Fe 


Equilibrium  Geometry3 

Lattice  Potential  Param°ters^ 

Defect0 

x(aQ) 

y(a0) 

z(a0) 

DLAT^kca^m0^  ^ 

“UT(ao''> 

none 

0.0 

0.0 

0.474 

127.9 

0.368 

step 

0.128 

0.0 

0.482 

111.0 

0.390 

kink 

0.128 

0.022 

0.482 

101.1 

0.408 

adatom-step 

0.128 

O.C 

0.482 

91.1 

0.429 

adatom 

0.0 

0.0 

0.474 

73.8 

0.473 

The  equilibrium  geometry  is  given  as  the  displacement  from  the 
unrelaxed  lattice  site.  The  z-axis  is  normal  to  the  surface, 
k  See  Equation  (1 ) . 

c  The  defects  are  discussed  in  Section  III. A. 
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TABLE  III.  LATTICE  POTENTIAL  PARAMETERS  FOR  FIVE 
SURFACES  OF  Ta. 


Defect0 

Equilibrium  Geometry3 

Lattice  Potential 

Parameters*5 

x(a0) 

y(a0) 

z<a0> 

DLAT(kcaVmol ) 

“LAT^o"'1 

none 

0.0 

0.0 

0.412 

245.7 

0.364 

step 

0.174 

0.0 

0.427 

214.9 

0.383 

kink 

0.132 

0.132 

0.430 

188.6 

0.401 

adatom-kink 

0.174 

0.0 

0.428 

172.6 

0.423 

adatom 

0.0 

0.0 

0.412 

140.1 

0.464 

3  The  equilibrium  geometry  is  given  as  the  displacement  from 

the 

unrelaxed  1 

at" ice  site.  The 

z-axis  is 

normal  to  the  surface. 

D  See  Equation  J ) . 

0  The  defects 

are  discussed  in 

Section  III. A. 

TABLE  IV. 

LATTICE 

POTENTIAL 

PARAMETER  FOR  FIVE 

SURFACES 

OF  Cr 

Equilibrium  Geometry3  Lattice  Potential  Parameters*3 


Defect0 

x(a0) 

yU0) 

z(a0) 

Dlat^03^/^0!  ) 

“latK  ^ 

none 

0.0 

0.0 

0.353 

117.3 

0.414 

step 

0.152 

0.0 

0.366 

102.6 

0.436 

kink 

0.115 

0.115 

0.368 

90.1 

0.461 

adatom-kink 

0.152 

0.0 

0.366 

82.5 

0.481 

adatom 

0.0 

0.0 

0.353 

67.0 

0.527 

The  equilibrium  geometry  is  given  as  the  displacement  from  the 
unrelaxed  lattice  site.  The  z-axis  is  normal  to  the  surface. 
See  Equation  (1 ) . 

The  defects  are  discussed  in  Section  1 1 1. A. 
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TABLE  V.  ATOM-ATOM  MORSE  PARAMETERS  FOR 
GAS  PHASE  DIATOMICS 


AB 

Morse  Parameters 

DAB(kcal/mol) 

“ab'V1' 

RAB<ao> 

H23 

109.0 

1.044 

1.401 

*2* 

228.4 

1.42 

2.07 

FeHb 

60.1 

0.74 

2.82 

a  Reference  12. 


Estimated  from  References  13  and  14. 


TABLE  VI.  ATOM-ATOM  MORSE  PARAMETERS  FOR  GAS-ATOM  . 

INTERACTIONS  WITH  SEMI-INFINITE  SURFACES3 


Metal 

Gas  Atom  B 

Morse 
Parameter 
.  Set 

Morse 

Parameters 

DgS(kcal/mol ) 

aBS(ao_1) 

RBS^ac^ 

Feb 

H 

4.50 

0.488 

6.07 

Fec 

N 

(1) 

6.30 

0.550 

6.75 

(2) 

7.55 

0.470 

6.40 

Tab 

H 

17.0 

0.841 

4.69 

Tac 

N 

(1) 

17.0 

0.841 

4.69 

(2) 

15.58 

0.580 

6.25 

Crc 

H 

(1) 

5.20 

0.495 

6.05 

(2) 

6.72 

0.564 

6.41 

Crc 

N 

(1) 

7.40 

0.545 

6.85 

(2) 

7.76 

0.462 

6.56 

a  For  those  cases  with  two  entries  the  values  presented  cover  a  range 


of  the  estimates  for  the  parameters. 
b  Reference  15. 
c  Estimated. 


TABLE  VII.  HYDROGEN  ADSORPTION  ON  SEMI-INFINITE 
SURFACES  OF  Fe . 


Defect^ 

Equilibrium  Geometry9 
x(aQ)  y(aQ)  z(aQ) 

Adsorption  Energy 
(kcaTOT) - 

Vibrational 

Frequency 

Normal  to  Surface 
(cm  ' ) 

none 

2.71  2.71 

3.60 

66.7 

690 

step 

3.11  2.71 

3.78 

61.0 

694 

kink 

3.05  3.05 

3.87 

56.4 

688 

a  The  origin  was  located  at  a 

lattice 

site  in  the  top  layer 

of  an  unrelaxed 

lattice.  The  calculations 

were  performed  with  the  top  layer  of  the  surface 

in  its 

relaxed  geometry. 

D  The  defects  are  discussed 

in  Section  III. A. 

TABLE  VIII. 

HYDROGEN 

SURFACES 

ADSORPTION  ON  SEMI-INFINITE 

OF  Ta. 

Defect^ 

Equilibrium  Geometry9 
x(aQ)  y(aQ)  z(aQ) 

Adsorption  Energy 
(kcal/mol ) 

Vibrational 

Frequency 

Normal  to  Surface 
(cm-1 ) 

none 

3.12  3.12 

1.39 

103.7 

1219 

step 

3.23  3.12 

1.45 

101.5 

1128 

kink 

3.20  3.20 

1.50 

99.8 

1075 

The  origin  was  located  at  a  lattice  site  in  the  top  layer  of  an  urirelaxed 
lattice.  The  calculations  were  performed  with  the  top  layer  of  the  surface 
in  its  relaxed  geometry, 
k  The  defects  are  discussed  in  Section  III. A. 


TABLE  IX.  HYDROGEN  ADSORPTION  ON  SEMI-INFINITE  SURFACES  OF  Cr 


Morse 

.  Parameter 

Equilibrium  Geometry3 

Adsorption 

Energy 

Vibrational 
Frequency 
Normal  to  Surface 

Defect0  Setc 

x(a0) 

yu0) 

z(a0) 

(kcal/mol ) 

(cm-T) 

none  (1) 

2.72 

2.72 

3.50 

74.0 

746 

(2) 

2.72 

2.72 

4.50 

75.0 

996 

step  (1) 

3.00 

2.72 

3.59 

68.4 

736 

(2) 

3.03 

2.72 

4.58 

69.5 

978 

defect  (1) 

2.99 

2.99 

3.67 

63.5 

728 

(2) 

3.02 

3.02 

4.65 

64.8 

961 

a  The  origin  was  located  at  a  lattice  site  in  the  top  layer  of  an  unrelaxed 
lattice.  The  calculations  were  performed  with  the  top  layer  of  the  surface 
in  its  relaxed  geometry, 
k  The  defects  are  discussed  in  Section  III. A. 

Q 

The  parameter  set  used  in  these  calculations  were  those  given  in  Table  VI. 
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TABLE  X.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE 
SURFACES  OF  Fe 


Defectb 

Morse 

Parameter 

Setc 

Equilibrium  Geometry3 

Adsorption 

Energy 

Vibrational 
Frequency 
Normal  to  Surface 

*U0) 

y(a0) 

z(a0) 

(kcal/mol ) 

~~~  (cm- ■') 

none 

(1) 

2.71 

2.71 

5.03 

79.7 

279 

(2) 

2.71 

2.71 

3.97 

119.7 

254 

step 

(1) 

3.14 

2.71 

5.10 

73.6 

273 

(2) 

3.11 

2.71 

4.07 

109.9 

251 

kink 

(1) 

3.12 

3.12 

5.17 

68.1 

257 

(2) 

3.09 

3.09 

4.16 

101.4 

248 

a  The  origin  was  located  at  a 

lattice 

site  in 

the  top  layer  of  an  unrelaxed 

lattice.  The  calculations  were  performed  with  the  top  layer  of  the  surface 
in  its  relaxed  geometry. 

13  The  defects  are  discussed  in  Section  III. A. 

c  The  par  ieter  set  used  in  these  calculations  were  those  given  in  Table  VI. 


38! 
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TABLE  XI.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE  SURFACES  OF  Ta. 


Defect*3 

Morse 

Parameter 

Setc 

Equilibrium  Geometry3 

Adsorption 

Energy 

Vibrational 
Frequency 
Normal  to  Surface 

x(a0) 

yu0) 

2U0) 

(koal/mol 

(cm-1 ) 

none 

(1) 

3.12 

3.12 

1.10 

140.1 

275 

(2) 

3.12 

3.12 

1.39 

103.7 

327 

step 

0) 

3.28 

3.12 

1.14 

134.1 

257 

(2) 

3.23 

3.12 

1.45 

101.5 

303 

kink 

(1) 

3.25 

3.25 

1.18 

129.0 

244 

(2) 

3.20 

3.20 

1.50 

99.8 

239 

a  The  origin  was  located  at  a  lattice  site  in  the  top  layer  of  an  unrelaxed 

lattice.  The  calculations  were  performed  with  the  top  layer  of  the  surface 

in  its  relaxed  geometry. 

b  The  defects  are  discussed  in  Section  III. A. 

£ 

The  parameter  set  used  in  these  calculations  were  those  given  in  Table  VI. 
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TABLE  XII.  NITROGEN  ADSORPTION  ON  SEMI-INFINITE  SURFACE  OF  Cr 


Defectb 

Morse 

Parameter 

Setc 

Equilibrium  Geometry3 

Adsorption 
, Energy 

Vibrational 
Frequency 
Normal  to  Surface 

x(a0) 

y(a0) 

Z(V 

(kcal/mol ) 

(cirri ) 

none 

(1) 

2.72 

2.72 

4.93 

89.9 

295 

(2) 

2.72 

2.72 

4.00 

129.9 

266 

step 

0) 

3.17 

2.72 

5.08 

82.8 

289 

(2) 

3.17 

2.72 

4.11 

119.0 

262 

kink 

(D 

3.15 

3.15 

5.15 

76.6 

283 

(2) 

3.14 

3.14 

4.20 

109.7 

259 

a 

The  origin  was  located  at  a  lattice  site  in  the  top  layer  of  an  unrelaxed 
lattice.  The  calculations  were  performed  with  the  top  layer  of  the  surface 
in  its  relaxed  geometry. 

13  The  defects  are  discussed  in  Section  III. A. 

c  The  parameter  set  used  in  these  calculations  were  those  given  in  Table  VI. 


t 
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TABLE  XIII.  Fe-Fe  MORSE  POTENTIAL  PARAMETERS  FOR 
A  FINITE  CLUSTER  OF  Fe  ATOMS 


Defect3 

Atom  in 
Cluster 

DggCkcal/mol ) 

aB$^ao  1 ^ 

*Wao^ 

none 

46 

10.00 

0.736 

5.34 

step 

38 

10.03 

0.735 

5.34 

kink 

33 

10.06 

0.735 

5.34 

adatom-step  34 

10.09 

0.734 

5.34 

adatom 

26 

10.14 

0.731 

5.35 

a  Defects  are  discussed  in  Section  1 1 1. A. 
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TABLE  XIV. 

ATOM-ATOM  MORSE  PARAMETERS 
ON  FINITE  Fe  CLUSTERS3 

FOR  HYDROGEN  ADSORPTION 

Defects5 

Atoms  in 
Cluster 

Morse  Parameters 

DgS(kcal/mol ) 

aBS^ao  ^ 

Wao) 

none 

30 

5.75 

0.478 

5.86 

step 

26 

5.74 

0.481 

5.91 

kink 

23 

5.77 

0.482 

5.91 

a  The  Morse  parameters  were  calculated  so  that  the  adsorption 
energies  and  vibrational  frequencies  in  Table  V.II  were 
reproduced  for  a  finite  cluster. 

5  The  defects  are  discussed  in  Section  III. A. 

TABLE  XV.  ATOM-ATOM  MORSE  PARAMETERS  FOR  HYDROGEN  ADSORPTION 
ON  FINITE  Ta  CLUSTERS3 


Morse  Parameters 


Defect5 

atoms  i n 
Cluster 

DBS(kcal/mol ) 

aBS^ao  ^ 

Wao^ 

none 

30 

0.841 

4.68 

step 

26 

17.12 

0.841 

4.68 

kink 

23 

17.11 

0.840 

4.  68 

a  The  Morse  parameters  were  calculated 

so  that  the  adsorption 

energies 

and  vibrational  frequencies 

in  Table  VIII 

were 

reproduced  for  a  finite  cluster. 

The  defects  are  discussed  in  Section  III. A. 
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TABLE  XVI.  ATOM-ATOM  MORSE  PARAMETER  FOR  HYDROGEN  ADSORPTION 
ON  FINITE  C r  CLUSTERS9 


Defects^ 

Atom  in 
Cluster 

Morse 

Parameter 

Setc 

Morse  Parameter 

DBS(kcal/mol ) 

aBS^ao  ) 

*Wao) 

none 

30 

(1) 

6.44 

0.486 

5.86 

(2) 

7.70 

0.559 

6.31 

step 

26 

(1) 

6.47 

0.484 

5.82 

(2) 

7.71 

‘  0.557 

6.28 

kink 

23 

(1) 

6.49 

0.483 

5.79 

(2) 

7.73 

0.555 

6.2 

d  The  Morse  parameters  were  calculated  so  that  the  adsorption 
energies  and  vibrational  frequencies  in  Table  XI  were 
reproduced  for  a  finite  cluster. 
b  The  defects  are  discussed  in  Section  III. A. 
c  The  parameters  sets  correspond  to  those  used  in  Table  VI . 
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TABLE  XVII.  ATOM-ATOM  MORSE  PARAMETERS  FOR  NITROGEN  ADSORPTION 
ON  FINITE  Fe  CLUSTERS* 


none  30  (1) 

(2) 

Step  26  (1) 

(2) 

kink  23  (1) 

(2) 


7.92 

0.545 

6.64 

9.79 

0.461 

6.16 

7.96 

0.544 

6.60 

9.87 

0.460 

6.12 

7.98 

0.544 

6.60 

9.93 

0.460 

6.11 

a  The  Morse  parameters  were  calculated  so  that  the  adsorption 
energies  and  vibrational  frequencies  in  Table  X  were 
reproduced  for  a  finite  cluster. 
b  The  defects  are  discussed  in  Section  III. A. 
c  The  parameters  sets  correspond  to  those  used  in  Table  VI. 
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TABLE  XVIII.  ATOM-ATOM  MORSE  PARAMETERS  FOR  NITROGEN  ADSORPTION 
ON  FINITE  Ta  CLUSTERS* 


Defectsb 

Atom  in 
Cluster 

Morse 

Parameter 

Sate 

Morse  Parameter 

DBS(kcal/mol ) 

aBS^ao  ^ 

RBS^ao^ 

none 

30 

(1) 

16.76 

0.579 

4.89 

(2) 

17.14 

0.841 

4.68 

step 

26 

(1) 

16.66 

0.577 

4.89 

(2) 

17.12 

0.841 

4.68 

kink 

23 

(1) 

16.59 

0.576 

4.90 

(2) 

17.11 

0.840 

4.68 

3 

The  Morse  parameters  were  calculated  so  that  the  adsorption 
energies  and  vibrational  frequencies  in  Table  XI  were 
reproduced  for  a  finite  cluster. 
b  The  defects  are  discussed  in  Section  III. A. 
c  The  parameters  sets  correspond  to  those  used  in  Table  VI. 


TABLE  XIX.  ATOM-ATOM  MORSE  PARAMETERS  IN  NITROGEN  ADSORPTION 
ON  FINITE  Cr  CLUSTERS3 


Defects^ 

Atom  in 
Cluster 

Morse 

Parameter 

Setc 

Morse  Parameter 

DBS(kcal/mol ) 

aBS^ao  ) 

RBS^ao^ 

none 

30 

(1) 

8.68 

0.539 

6.73 

(2) 

10.26 

0.453 

6.30 

step 

26 

(1) 

0.72 

0.538 

6.70 

(2) 

10.31 

0.451 

6.26 

kink 

23 

(1) 

8.72 

0.536 

6.68 

(2) 

10.39 

0.451 

6.239 

'f  •  a  The  Morse  parameters  were  calculated  so  that  the  adsorption 

energies  and  vibrational  frequencies  in  Table  XII  were 
reproduced  for  a  finite  cluster, 
k  The  defects  are  discussed  in  Section  III. A. 
c  The  parameters  sets  correspond  to  those  used  in  Table  VI. 


Defects  are  discussed  in  Section  III. A. 
Exoergicity  for  the  Rideal  process. 


RAJECTORY  RESULTS  FOR  ATOM  RECOMBINATION  ON 
ETAL  SURFACES  BY  A  RADEAL  MECHANISM3 
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TABLE  XXII.  RESULTS  OF  TRAJECTORY  CALCULATIONS  FOR  NITROGEN  ATOM 
RECOMBINATION  ON  Fe  SURFACES  FOR  NON-NORMAL  INCIDENCE 
OF  THE  GAS  ATOM3 


Defect 

et 

(kcal/mol ) 

6 

♦ 

Pr 

AE$ 

(kcal/mol) 

aes/aer 

none 

9.24 

O 

o 

CO 

o 

O 

co 

.63 

46 

.31 

9.24 

30° 

150° 

.35 

57 

.39 

step 

4.62 

o 

o 

CO 

00 

o 

o 

.89 

45 

.29 

kink 

4.62 

30° 

CO 

o 

o 

.92 

51 

.32 

a  All  trajectories  are  for  Morse  parameter  set  (1)  and  a 
surface  temperature  of  1500K. 
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Previous  models  for  gas-- surface  Interactions  have  utilized  either  pairwise  forces  or  semi-empirical  potentials  formulated 
in  terms  of  perfect  two-dimensional  su. faces.  Here,  we  combine  the  pairwise  model  with  existing  semi-empirical  models, 
retaining  the  realistic  chemical  forces  of  the  semi-empirical  model,  but  having  direct  applicability  to  surfaces  with  non¬ 
periodic  structures, 


1.  Introduction 

In  a  previous  paper  [1],  we  presented  a  model  po¬ 
tential  that  has  the  appropriate  characteristics  for  de¬ 
scribing  the  reactions  of  diatomic  molecules  with  solid 
surfaces.  One  of  the  primary  objectives  of  this  line  of 
research  has  been  to  model  in  a  realistic  way  the  dyn¬ 
amics  of  chemical  reactions  occurring  on  and  with  solid 
surfaces.  For  this  end,  a  model  potential  must  have 
three  properties:  (1)  It  must  be  realistic,  reproducing 
the  experimental  binding  energies,  activation  barriers, 
etc.,  so  far  as  they  are  known;  (2)  it  must  be  flexible, 
capable  of  modelling  a  variety  of  different  potentials 
by  the  same  overall  procedure;  and  (3)  it  must  be  trac¬ 
table,  such  that  forces  experienced  by  each  atom  can 
be  rapidly  computed  and  classical  trajectory  studies  can 
be  performed.  The  London-Eyring-  Poianyi-Sato 
(LEPS)  [2 j  seemed  to  meet  these  three  criteria  and 
has  been  used  as  the  starting  point  for  a  variety  of 
studies  of  heterogeneous  reaction  dynamics  [3  j . 

While  the  original  LEPS  model  treated  the  solid  surface 
as  rigid  [2 j ,  it  was  soon  generalized  to  permit  oscilla¬ 
tions  of  the  surface  atoms  [4  j ,  generalized  Langcvin 
oscillutors  [5],  and  reactions  with  the  solid  itself  [!]. 
However,  the  model  potential  as  implemented  thus 
far  has  only  been  used  for  surfaces  with  perfect  two- 
dimensional  symmetry  (apart  from  local  oscillations 

*  This  research  was  sponsored  by  the  U.S.  Army  Research 
Office  under  Contract  Number  DAAG29-77-0047. 


or  displacement  from  equilibrium  lattice  positions), 

In  this  note,  we  point  out  that  the  LEPS  potential  can 
be  re-formulated  in  such  a  way  that  imperfect  crystal? 
can  be  modelled  keeping  the  same  realism,  flexibility, 
and  tractability  that  has  been  so  useful  in  these  recent 
studies.  There  is  considerable  evidence  that  many  cau- 
lystic  reactions  proceed  on  steps  [6]  and  amorphous 
metals  have  many  interesting  surface  chemical  proper¬ 
ties  [7].  Therefore,  It  is  important  to  add  ways  to 
study  such  systems  to  our  arsenal  of  computational 
procedures. 

In  this  note,  wc  alter  the  modified  LEPS  potential 
function  to  allow  the  treatment  of  non-rigid  surfaces 
with  arbitrary  geometries.  This  could  include  genuine!;- 
amorphous  materials,  polycrystallino  samples  with  step 
kinks,  etc.  The  alteration  Is  straightforward.  Previously 
[2]  the  interaction  of  one  gas  atom  with  the  solid 
surface  was  assumed  to  be  a  Morse  function  in  the 
direction  perpendicular  to  the  surface.  The  Morse 
parameters  were  allowed  to  be  periodic  functions  of 
the  X  and  Y  distances  over  the  surface.  The  values  of 
the  parameters  were  fixed  to  those  values  which  gave 
agreement  with  known  experimental  or  theoretical 
values  for  binding  energies  of  the  fragment  to  the  sur¬ 
face.  We  call  this  potential  i3pg.ri0dic,A  tinderstanding 
that  it  describes  the  interaction  of  a  single  atom  (A) 
with  a  solid  smface,  periodic  in  X,  Y.  The  final  LEI’S 
model  was  given  as  a  functional  of  Cpor iodic, A  ai,L*  l'ie 
potential  describing  (he  interaction  of  the  two  gas 
atoms  jj.  Therefore,  the  LEPS  formulation  for  the 
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periodic  surface  is 

'  pet iodic  n  periodic, A  '  Qpcriodic.il '  ^A.B^  ' 

Correction  terms  have  been  added  unto  tltis  basic  mod¬ 
el  potential  to  allow  for  motion  of  the  surface  atoms, 
erosion,  etc.,  but  they  are  not  essential  to  our  present 
discussion. 

The  present  modification  replaces  (Aperiodic  a 
sum  of  pairwise  forces  between  the  gas  atoms  and  each 
surface  atom,  but  retains  the  lunctional  dependence 
lii  cq,  (i ).  That  is,  if  S  denotes  the  set  of  atoms  in  the 
solid,  the  usual  pairwise  model  for  atom-surface  inter¬ 
actions  is  simply 


Qpairwisn.A  B  ^/A  >  (2) 

whore  q!A  is  the  two-body  force  betsveen  solid  atom  i 
and  gas  atom  A.  A  fully  pairwise  description  of  the 
interaction  between  diatomic  molecule  AI3  and  the 
solid  is  then 

^pairwise  ~  QpnlrwUe,A  +  Qpalrwise.B  +  ^A,B  *  (3) 

As  discussed  elsewhere  (1],  tius  is  r.ot  expected  to  be 
satisfactory  for  many  cases  of  interest.  However,  the 
modification  proposed  here  is  simply  to  use  ^pairwise 
to  0 )  in  place  of  <2p<rrJodio  8ivinS 

^amorphous  “  ^^Qpaitwisc, A  •  Qpalrwisc.B  •  ^A.B^  •  W 
Thus  the  entire  LEPS  formulation  given  previously 
follows  through  changing  simply  the  form  of  the  under¬ 
lying  component  Interactions,  The  form  of  the  potential 
for  diatomic  AB  on  surface  5  is 


^LEPS  *  t/AB+  ^AS*  ^BS  ^ 

“  ^AS  +  ^US)2  "  A  As('4AS  +  /1Bs)l 

where 

Da  r> 

^AB  1  4(TTa^)  {(3  '  AAb)«P  l-2aABG’AB-/?AB)J 

-(2+6AAg)exp  1"0ab(^AB  ~  /?ab)U  ’  (6) 

•^AB 

■^AB  -4075^)K,  +  3^AIj)cxP  [“‘“aU^AB'^Ab)} 
-(6+2AAU)cxp  l-aAB(rAfl  -  /*AB)]}  ,  (7) 


V 


AS 


I 


A  i 


A(k 


X{(d  +  ^AS)exp  I~2aA/(rA/  - RAi)) 

-  (2  +  6Aas)  exp  [-«A/(rA/  -  RAi)}}  , 

Xdai 


(8) 


‘■S 4(1 +  .!*)& "A/ 

X  ((1+  2AAS)exp  (-2aA/(rA/-flA/)J 
-  (6+  2AAS)c.xp  [-orA<(rA/ -  jRa/)]}  , 


(9) 


and  2)pQ,  erpQ,  and  R pq  are  the  dissociation  energy, 
Morse  parameter  and  equilibrium  distance  for  the  tsvo- 
body  interaction  between  P  and  Q.  $  and  ,-lBS  are 
derived  from  eqs.  (8)  and  (9)  by  substituting  B 
for  A.  Apq  is  the  Sato  parameter  used  to  generate  a 
variety  of  surfaces.  Only  a  change  in  the  expressions 
for  the  two-body  atom-surface  Coulomb  and  exchange 
integrals  has  been  required.  The  use  of  the  pairwise 
additive  atom -solid  potential  as  a  single  two-body 
potential  is  consistent  with  the  requirement  (imposed 
during  the  derivation  of  the  modified  LEPS  potential) 
(2]  that  the  surface  acts  as  an  extended  source  of  one 
and  two  electrons,  as  required,  for  bonding  the  incident 
fragments. 

In  fig.  1  we  compare  the  periodic  potential  from  ref. 
[1  ]  (fig.  3b)  for  Il2  interacting  with  the  (00’)  face  of 
iron  with  Kumcirpllous  obtained  from  conside  ing  the 
interaction  with  a  cluster  of  33  Fe  atoms  (fig.  2).  The 
pairwise  forces  qiA  in  eq.  (2)  have  the  form  of  Morse 
potentials 


qlA  »  D o  exp  [-(m/r0)  (r)A  -  r„)] 

X  {exp  [~(m/r0)(r<A  -  r0)]  -  2} ,  (10) 

where D0  =  4.6  kcal,  r0  =  3.1 1  A.  m  =  2.96.  We  see  in 
fig.  1  that  the  potentials  have  some  differences  quantita¬ 
tively,  but  have  the  sane  qualitative  structure.  That  is, 
there  is  an  activation  barrier  lo  adsorption  approximately 
in  agreement  with  experiment  [1 ,8]  and  no  molecular 
adsorption  is  present.  For  large  H-H  separations,  when 
both  atoms  lie  on  the  surface,  there  are  some  effects 
from  the  finite  size  of  the  cluster  and  Hie  neglect  of 
continuum  contributions  [9],  Also,  the  H?  -Fc  inter¬ 
action  more  rapidly  decreases  for  the  cluster  than  for 
the  infinite  surface  as  ’i2  recedes  from  the  surface. 
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Fig.  1.  F.quipotential  contours  (in  «V)  fat  the  approach  of  Hj 
towards  the  Iron  surface.  The  H-H  bond  is  puralicl  to  the  sur¬ 
face,  the  bond  midpoint  is  centered  over  a  ICN  site,  and  oriented 
in  the  direction  of  adjoining  5CN  sites  (defined  In  fig.  2).  This 
path  lends  to  the  smallest  barriers  to  adsorption  as  well  a'-  the 
most  tightly  bound  adsorbed  atoms,  (a)  Aperiodic  ltorn  i!F 
(1);  (b)  •''amorphous  from  cq.  (4). 


However,  at  present  there  is  no  way  to  determine  from 
experiment  which  fit  is  preferable,  so  we  forgo  exten¬ 
sive  variation  of  parameters  to  produce  a  fully  adjusted 
surface. 


Fig.  3.  Interaction  generated  from  •■'pairwise  using  On  The 
othc,  conditions  are  lire  same  as  fig.  1. 
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In  conclusion,  the  substitution  of  a  pairwise  addi¬ 
tive  fragment-surface  interaction  potential  for  the 
P<cviously-usod  periodic  potential  seems  to  yield  a 
use t'ul  potential  function  for  the  study  of  gas- 
jimuphous  surface  interactions.  It  retains  most  of  the 
flexibility  of  the  previous  model  while  providing  cor¬ 
rect  functional  forms  for  several  new  cases  of  interest 
including  amorphous  solid  structures,  energy  transfer 
to  the  surface,  and  erosion. 
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